BARI-TH/97-266 
May 1997 



Violations of Local Duality 
in the Heavy Quark Sector 

P. Colangelo a , C.A. Dominguez 6 , G. Nardulli a ' c 

a Istituto Nazionale di Fisica Nucleare, Sezione di Bari, Italy 
Institute of Theoretical Physics and Astrophysics, University of Cape Town, 

Rondebosch 7700, South Africa 
c Dipartimento di Fisica, Universita di Bari, Italy 



Abstract 

We examine the origin of possible failures of the quark-hadron local duality. In 
particular, we consider a correlator of two currents comprising heavy quark oper- 
ators, and compare the evaluation by the Operator Product Expansion with the 
result obtained inserting an infinite set of hadronic states, in the iV c —* oo limit. 
Whereas the smeared spectral functions agree with each other, their local behaviour 
is different. The difference manifests itself in the form of a term (D(A/e) (e being 
the residual energy) that is not present in the expansion in powers of A/e predicted 
by the OPE. 



1 Introduction 



Quark-gluon/hadron duality represents a key concept in the theoretical description of 
inclusive hadronic processes. By its use, high energy processes such as, e.g., e + e~ an- 
nihilation, can be computed in terms of hadronic matrix elements of operators in an 
expansion, the Operator Product Expansion (OPE), whose leading term is represented 
by the perturbative Quantum-Chromo-Dynamics expression. As a matter of fact, duality 
provides the tool to extrapolate from the deep Euclidean region, where the OPE is defined, 
to the Minkowski domain, where physical observables are measured. The basic assump- 
tion underlying this trade-off between hadronic quantities and quark-gluon amplitudes is 
that, at high energy and/or momentum transfer, the hadronic behaviour should be well 
described by quark and gluon interactions, provided that sufficiently inclusive variables 
are considered. 

In the last few years this approach has been extended to the decays of hadrons con- 
taining one heavy quark, exploiting the presence in these systems of a large parameter, 
the heavy quark mass uiq, that can be taken, at least formally, infinitely large: niq — > oo 
0. In this extension, however, a clear distinction must be maintained between semilep- 
tonic and nonleptonic decays. As a matter of fact, in general, the hadronic and the 
OPE amplitudes cannot be identical even at very high momentum transfer, due to the 
different structure of their singularities: multiparticle thresholds in the former case and 
quark-gluon production thresholds in the latter. For semileptonic heavy hadron decays, 
however, this difficulty does not prevent duality to hold; as a matter of fact, in computing 
the semileptonic inclusive width, one has to integrate over lepton variables, which amounts 
to a smearing of the OPE width. The equality between smeared OPE and physical, i.e. 
hadronic, widths is sometimes referred to as global duality (assumed to be distinct from 
local duality, i.e. without smearing) and it is generally believed that global duality holds 
between quark-gluon and hadronic cross sections 0; for B and decays it has been 
proved explicitly in [[| to two orders in the 1/toq expansion and to the first order in a s , 
in the so-called small-velocity (SV) limit {m h ,m c » m h — m c » Kqcd) [§]• 

For nonleptonic heavy hadron decays, on the other hand, there are no lepton momenta 
to be integrated and, therefore, one cannot use global duality to prove the identification 
of OPE and hadronic observables. Therefore local duality must be assumed, which is a 
stronger hypothesis. For this reason, the validity of OPE for the computation of nonlep- 
tonic heavy hadron inclusive decays appears more debatable. For example in || it has 
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been argued that the discrepancy between the OPE prediction > 0.9 Q and the 

experimental result — — b — = 0.78 ± 0.04 || might be solved assuming a violation of local 

^~ v d) 

duality with the appearance of a O(^) correction not predicted by OPE. In this context, 
it may be useful to observe that in |10[ local duality has been proved for nonleptonic Af, 
and B inclusive decays, for the first two terms in the 1/rrib expansion, and at the order a s 
in the perturbative expansion, assuming factorization of the weak amplitudes. Also this 
result, however, holds in the SV limit, i.e. in a kinematical regime which is far off from the 
physical one; therefore, in absence of a general proof, the failure of local duality, with the 
appearance of O(^) corrections to OPE, remains a logical possibility and constitutes a 
simple explanation of the A?, lifetime data. 

It would be clearly extremely important to understand, in this scenario, the origin 
of this possible ^ correction to the Operator Product Expansion, but this appears, at 
the moment, a formidable task since, at the present stage of (analytical) understanding of 
non-perturbative QCD, no first-principle answer can be given to the question of how good 
duality really is and which violations it may suffer. A more modest, but, nevertheless, 
potentially instructive aim could be to prove the existence of a violation of local duality 
and the emergence of a ^ correction in some definite model. This is the main purpose 
of the present letter. 

There have been some other recent studies on the validity of the quark- gluon/hadron 
duality in connection with non-perturbative QCD applications in the heavy quark sector 
PP , [T2|l . Also these works are motivated by the hope to improve the precision with which 
various observables are determined using the OPE coupled with duality. For example, 



the model studied in [1 1J is based on an Ansatz for the evaluation of the two-point func- 
tion involving the difference between scalar and pseudoscalar heavy-light currents. In 
the chiral limit (for the light flavour) and the infinite mass limit (for the heavy flavour), 
the coefficients of the OPE for this correlator can be calculated analytically. A spectral 
function model (in the time-like region) is then expanded in a power series (in the space- 
like region) and compared with the exact result. From this comparison one can gauge 
the validity of (local) duality. In this letter we consider the same two-point function and 
we analyze two models: the first model, to be studied in the next Section, represents a 



slight modification of ]IT[ and leads to a maximal violation of duality. The second model, 
discussed in Section 3, uses results for the hadronic evaluation of the two-point function 
1 For a discussion see ||, 0, |). 



obtained by a constituent quark approach; therefore it represents an improvement as com- 
pared to the previous Ansatze that are ad hoc assumptions without physical justification. 
This calculation explicitly shows a violation of local duality in the form of an unexpected 
O(-) correction to OPE. 

2 Mathematical models for the two-point function 

We begin by considering the following two-point function: 

1%) = l -f dxe*** < 0|T(J 5 (x)4(0)) - T(J p (x)JU0))\0 >= U s (q) - U P (q) (1) 
where the scalar (J s ) and the pseudoscalar (Jp) currents are 

Js(x) = Q(x)q(x) (2) 

J p [x) = Q(x)i^ 5 q(x) (3) 

and Q(x), q(x) are heavy and light quark operators, respectively. This two-point function 
is particularly simple, as in the chiral limit m q — > it vanishes in perturbation theory. 
Hence, 11(g) is entirely non-perturbative. Additional simplifications take place in the 
infinite mass limit tuq — > oo, where it is convenient to write 

q» = (m Q -e,q = 0) . (4) 
In this limit the correlator becomes a function of e, i.e. 

H(c) = \l +0 °d^~ tT ®(r) (e>0) (5) 



where $(r) is sometimes called the non-local quark condensate |T3 

$(v^2) = <0\q( x )e l9s fo d y* A ^q(0)\0> . (6) 

In the limit e >> Aqcd, the OPE expression for 11(e) is given by 

< Q\qq\Q > 



UopE{e) 



4e 



2 4 fi 

mi rrin rrvf. 
8e z e 4 e b 



where < 0|gg|0 >= (-240 MeVf, m 2 = < °l^^^|0 > = .8 ± 0.2 GeV 2 . The 

mi K ' ' < 0|gg|0 > 

positive coefficients Cm above, depend on the actual form of the non local condensate. 
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(7) 



We note explicitly the alternating signs in Eq.(|7|) and the absence of even powers of e 1 ; 
both these features are consequences of the general principles on which the OPE is based 



lTJ. In particular, the absence of a D = 4 term, proportional to the gluon condensate 
< a s G 2 >, is due to the limit mq — > oo, m q — > 0. 



Let us now briefly review the model proposed in []TTj] to analyze duality violations. It 
is given by 

< 0|gg|0 > „,e + A 



n(c) 



where A is a parameter and 



i3(z) 



if, 



4A 



z + 1 



-/?(- 



2A 



(9) 



ip(z) being the logarithmic derivative of the Gamma function. The model correlator 
admits the following series expansion 



rn 



< 0\qq\0 > 
2A 



E 



" G Ae + 2j + l 



(10) 



where A = 1/A. From this expression for 11(e) we get, by the Mellin transform, 

< 0|gg|0 > 



$1 



T 



cosh(Ar) 

For e — ► oo, 11(e) has the asymptotic expansion 

< 0|gg|0 > 



n e 



4e 



E E 



A 



2 a 



2n" 



where E 2n are Euler numbers (Eq = 1, E 2 



n=0 

-1, E4 = 5, Eq 



:2n 



(11) 



(12) 
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Comparison of this result with that of the OPE, Eq.(|7|), indicates that this model 
is able to reproduce the right power structure of 11(e). The associated spectral density 
corresponds to an infinite number of equally spaced poles located along the negative e 
axis, these poles having residues alternating in sign. After smearing, the spectral function 



looks rather reasonable |TTJ. However it is easy to show that this result strongly depends 
on the Ansatz @. As a matter of fact, let us consider the slightly modified model 

-IV 



11(e) = const, x 



Pi Ae + 2j 



(13) 



instead of (10). This implies 



$(t) = const, x [tanh(-) — 1] 

A 



(14) 



In the limit e 



oo, this 11(e) becomes 

n (e) e^oo const - x 



1 + E 



G 



n=0 



2n+l 



(15) 



A comparison with the OPE result, Eq.(^), shows that except for the first term, the power 
structure is wrong, so that one obtains a maximal violation of duality. This model is still 
quite realistic. The spectral function is also made of an infinite number of zero-width 
resonances, with alternating sign residues, and located along the negative e axis. The 



only difference with the model of |TlJ is a shift in the location of the poles, due to the 
absence of the unit factor in the denominator of Eq . (|13|) . 

From this example we learn that the agreement found in [TT[| between the hadronic 
evaluation of the correlator and the OPE result might be fortuitous; in order to get a 
deeper understanding we compute now the correlator Eq.(|l]) using a more realistic model 
for the spectral density. 



3 A realistic model for the hadronic correlator 

We evaluate the two-point function ([!]) by inserting hadronic states between the currents. 
In general this can only be done in some approximation: we choose to insert an infinite 
number of states, but in the iV c — > oo limit (N c = number of colours), where the surviving 
contributions are J p = + and 0~ single particle states, contributing respectively to the 
scalar (lis) and pseudoscalar (lip) part of II. By denoting \S n > and \P n > these states, 
with masses M$ n and Mp n respectively, we define the current-particle matrix-elements: 

Ml 

<0\J s \S n > = -^f Sn (16) 

m Q 

M 2 

< 0\J P \P n > = — ^/ Pb . (17) 

m Q 

In the mQ — > oo limit we have 

M Sn = m Q + 5 Sn + 0(^) (18) 

M Pn = mQ + 5 Pn + 0( — ) . (19) 

m Q 

The binding energies 6s n and 5p n can be obtained by solving the wave equation 



-V 2 + m 2 Q + .^-V 2 + m 2 + V(f)\ * n (f) = M n ^ n (r) (20) 
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in the limit itlq — > oo,m q — > 0. Assuming a central potential V(f) = V(r) we can write 

. The scalar particles are obtained for i = 1 (P-waves) and the 



Y fm (f, 



u< n) (r) 



pseudoscalar particles for £ = (5— waves). The corresponding equation [15j for <55 n , <5p n 
is given, in this limit, by 



V(r) - 6, 



(n) 



i \ 2 f°° f°° i \ 

uT ''Or) H — / c?r' / dkkxe(kr)xe(kr')u" Or') 

71 JO JO 







(21) 



where the relation between 5s n , 5p n in (p~8| , pT9|) and 5^ is as follows: 5^ n) = 8 Pn , 5{ n) = 5 Sn ; 
moreover, Xl{ x ) = x Ji( x ) {je are the spherical Bessel functions). 

V{r) is the static interquark potential; several forms have been studied in the literature 
which reproduce the experimental spectrum of the heavy Qq mesons. In general they 
behave linearly V(r) ~ fi 2 r for r — > oo and have a coulombic behaviour (V(r) ~ — ) at 
small distances. In order to simplify our discussion we assume 



V(r) 



fi 2 r 



(22) 



with constant fi (string tension). This form is adequate for our aims because the coulombic 
part of the potential mainly affects the first states (n = 0) and becomes negligible with 
increasing n. 0. 

In order to solve Eq.(21) for any integer n (n = 0,1...), we work in the WKB ap- 
proximation [16]. The spectrum and the wavefunctions are obtained by the usual WKB 
procedure O, |16| ; the spectrum is given by 



f(n) 



fiJn(2n + e + -) ; 



(23) 



the wave function, for r < r = reads as follows: 



u 



f\r) = vV^Q Xe 



(r < r ) , 



(24) 



whereas for larger values of r it decreases exponentially. Let us observe explicitly that 
the factor JrriQ arises from the covariant normalization condition 



dr \tyjr)\ 2 = 2M„ 



(25) 



2 For example, using the Richardson potential and solving numerically the wave equation |iq| one 
obtains a coulombic correction to the meson mass as follows: AM n /M n = 7.9%, 5.6% and 4.4% for the 
radial quantum number n — 0, 1, 2 respectively. 
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Let us now turn to the coupling fs n = fi and fp n = f^ 1 ' . From the expression [117 



r(n) 



(n) _ /3 1 



2 7rikf n ^ jo 



with the function (A;) defined as 



c "1 / 
Ut (k) = dr xe{kr)ue(r) , 
Jo 



one obtains 



(<<) 



\ 



(26) 



(27) 



(2? 



7T M n> £ 

Let us now reconsider the correlator Eq. ([!]). In the iV c — > oo limit the products of currents 
in Eq.(l) are saturated by single particle states with J p = + and 0~, respectively. 
Therefore one can write: 



-i oo 

E 



5m 



n=0 



[ft 



Jo 



in) 



i v( n ) i r( ra ) 



(29) 



which, in the limit tuq — > oo, becomes 



,2 oo 



3/r _ 
8^ to 



(n) 



v(«) 



- (-l)" + Vn + 3/2 



o e + n Jiz(n + 3/2) 



(30) 



Before continuing in the analysis, a comment is in order on the model. The choice of 
considering only single particle states between the currents in the correlator (i.e. taking 
the N c — > oo limit) is not too restrictive, as it can be shown by considering, for example, 
the imaginary part of the correlator of scalar currents IT 5 in Eq.(jl]) (the same result is 
obtained taking Hp). By computing the imaginary part of the quark loop diagram one 
obtains, for E = — e >— > +00 and at the leading order in E: 

3E 2 



The resonance model gives 



ImU 



had 



ImH° s PE {E) 



im- - 



8tt 



n -E + 5^ - te 



|E[A ( 1 M^E-Sl 



3^" 



n=0 



E — fiJ2ir(n + 7/4) 



(31) 



(32) 



Im Tlf d (E) looks very different from Im Ug PE (E) in Eq.([H]), due in particular to the 
presence of the infinite set of Dirac 5 functions in Eq.(32). On the other hand it is well 
known that a comparison should be made only after a smearing of these distributions @. 
A convenient way to do this in the present case is make the approximation 



E 



dn 



(33) 



Hence, we obtain 



Im U h s ad 



3fj?E 



3E 2 



dn5 



E 



8tt 



ImTi° s PE {E) . 



(34) 



Therefore Hs(E) and Up(E) satisfy global duality at least at the leading order for E — > oo. 

In the sequel we shall discuss the next-to-leading contributions that are responsible 
for the difference II = lis — lip; for the time being let us comment on the value of the 
parameter /i and the approximate choice Eq.(22) for the potential. Assuming that the 
WKB approximation works reasonably well already for the the first quantum number 
n = 0, we obtain 



/' - \l^( m B- m b) ; 



(35) 



for 777,6 — 4.6 — 4.7 GeV this gives /i ~ 300 MeV. The first (n = 0) meson state may be 
sensitive to the small distance, coulombic part of the interquark potential that we have 
omitted in Eq.(22); however, as far as we neglect consistently 0(a s ) terms also in the 
OPE counterpart of II , this approximation is reasonable and the duality should be 
equally valid. 

Let us now apply the same resonance model to the correlator 11(g) = Ils(g) — Up(q) 
in Eq.([IJ). When we compute the correlator inserting hadronic states, from Eq.(|3"0D we 
find (£ = -e > 0): 



R had (E) = ImW ad (E) 



3fi 2 E 



£ (-!)"<$ 



n=0 



E — ir(n 



(36) 



On the other hand, from the first four terms in Eq.(^), we have: 



R 



OPE 



(E) 



ImU 

7i 



OPE 



(E) 



< 0\qq\0 > 



5(E) - ^5"(E) + c^SW(E) - c^5^\E)... 
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5 61 

C4 and cq are unknowns: in the model of [TJ] they are given by C4 = — , cq = . Also in 



64' 512 

the present case the two expressions look very different, but they can be compared after 
an appropriate smearing; following we consider 

R had (E, A) = — dE' — ^— $P A9 (38) 



7T J (E — E' f + A 2 

to a) = f ( /_°;y;» A2 . m 

In the limit A — > 0, -R(-E, A) — >■ R(E); however one has to impose A >> Aq C d @- As a 
matter of fact, 

R(E, A) = — [U(E - iA) - IL(E + iA)] (40) 

and, for A >> Aqcd, we are far enough from the singularities for OPE to be valid and 
duality to hold. 

For very large E the two expressions should coincide regardless of A. From previous 
equations one finds that this request implies the relation 

, , 0.51 u 3 , , 

< |gg|0 >= -3— -5=, (41) 

2 ^JTX 

and, for < 0|gg|0 >= (-240MeV) 3 , this implies fi = 317 MeV; it can be noted that 
this number agrees nicely with the previous result Eq . (|35|) . For smaller values of E, i.e. 
1 < E < 20 GeV, the numerical results concerning the comparison between fl38|) and (|39D 
are reported in Fig. 1, where we plot the ratio: 

jDhad 

P =BOPE ( 42 ) 

of the smeared quantities as a function of E for several values of A. We observe that, as 
expected, the agreement between R had and R OPE increases with increasing A; in particular 
for A = 3.0 GeV the difference does not exceed 20%. In any case, for any value of the 
smearing parameter A, the ratio (^) tends to unity for large E values. 

A further test of the soundness of our approximation is provided by the evaluation of 

F = ^T b f B (43) 



in the limit m& — > 00, using the WKB formula (|28|) with n = I = 0. We get F = 0.26 
GeV 3 / 2 to be compared with the QCD sum rule result, obtained without a s corrections, 
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i.e. in the same approximation used in this work: F = 0.30 ± 0.05 GeV 3 / 2 ||TB[ (including 
a 9 corrections, one would get F = 0.41 ± 0.04 GeV 3 / 2 ). 

In conclusion we have tested that the resonance model we have proposed satisfies 
global duality: the smeared imaginary parts of the correlator, when computed by OPE 
or by hadron states agree with each other. When we consider local duality, however, 
the expressions are significantly different. As a matter of fact, starting from Eq. (|3"0"|) , we 
obtain, by the Mellin transform: 



with 



nfr) 



3/i 



IL(t) 



3 oo 



1 

2vri 



u+ioo 



de e eT U(e) 



^(-l) n+1 v /vr(n + 3/2) e -W-(™+ 3 / 2 ) 



(44) 
(45) 



n=0 



IT can be expanded for small r; we find numerically 



U(r) 



-0.51 + 



0.46 



/i 2 7rr 2 



Since 



we find 



n(e) 



dre- Te U(T) 



U had (e) 



0.51/i 3 l 
8^ 7 



1-^ + 0.93^?... 

e e 2 



ml 



(46) 



(47) 



(4£ 



The factor multiplying - coincides numerically with <0 ^ 4 g ^ 0> for [i = 317 MeV (see eq . (|4lD ) . 
As for the mass parameter mo, numerically we find mo = 560 MeV. In conclusion we get 
(6<0): 



n 



had i 



< 0|gg|0 > 
4e 



m 



T-,2 



ffl, 



+ 0.93^ + 



(49) 



A comparison between Eq.(^) and Eq . (|49|) shows a violation of local duality which mani- 
fests itself in the form of an unexpected term in the - asymptotic expansion. 

e 



4 Conclusions 



Studying simple correlators of quark currents, we have found, in a particular quark model 
and in the iV c — > oo limit, an explicit example of violation of local duality which occurs in 
a condition where global duality is verified. The difference between the hadronic and the 
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OPE spectral functions is made evident by the A/e term which is absent in the expansion 
predicted by OPE. 

The calculation in QCD of correlators of quark operators between hadronic states, as 
required for the evaluation of, e.g., the B d and A b inclusive decay widths, is, of course, 
beyond the aims of the present analysis. However, the example described in this paper 
supports the suggestion that the Af, lifetime data could find an explanation in the presence 
of a 1/rriQ correction not included in the usual QCD-OPE expansion. 



Acknowledgments 

We thank N. Paver for useful discussions. The work of (CAD) has been supported in part 
by the FRD (South Africa). 



11 



References 



[I] J.Chay, H.Georgi and B.Grinstein, Phys. Lett. B 247 (1990) 399; I.Bigi, M.Shifman, 
N.Uraltsev and A.I.Vainshtein, Phys. Rev. Lett. 71 (1993) 496; B.Blok, L.Koyrakh, 
M.Shifman and A.I.Vainshtein, Phys. Rev. D 49 (1994) 3356; A.V.Manohar and 
M.B.Wise, Phys. Rev. D 49 (1994) 1310; A.F.Falk, M.Luke and M.J. Savage, Phys. 
Rev. D 523 (1996) 2491. 

[2] E.C.Poggio, H.R.Quinn and S.Weinberg, Phys. Rev. D13 (1976) 1958. 

[3] C.G.Boyd, B.Grinstein and A.V.Manohar, Phys. Rev. D 54 (1996) 2081. 

[4] M.B.Voloshin and M.Shifman, Yad. Fiz. 47 (1988) 801 (Sov. Journ. Nucl. Phys. 47 
(1988) 511). 

[5] G.Altarelli, G.Martinelli, S.Petrarca and F.Rapuano, Phys. Lett. B 382 (1996) 409. 
[6] M.Neubert and C.T.Sachrajda, Nucl. Phys. B 483 (1997) 339. 
[7] PColangelo and F.De Fazio, Phys. Lett. B 387 (1996) 371. 

[8] I. I.Bigi et al., in "B decays", S.Stone ed., 2nd edition, World Scientific (Singapore) 
(1994) 132. 

[9] J.Richman, talk given at the 28th International Conference on High Energy Physics, 
ICHEP96, Warsaw, Poland, 25-31 July 1996, to appear in the Proceedings. 

[10] A.Deandrea, N.Di Bartolomeo, R.Gatto and G.Nardulli, Phys. Lett. B 398 (1997) 
353. 

[II] M.Shifman, Theory of Preasymptotic Effects in Weak Inclusive Decays, in Proc. of 
the Workshop Continuous Advances in QCD, ed. A.Smilga (World Scientific, Singa- 
pore, 1994), p. 249; M.Shifman, in Proc. of the Johns Hopkins and Pascos Workshop: 
Particle, Strings and Cosmology, ed. J. Bagger (World Scientific, Singapore, 1996). 

[12] B.Chibisov, R.David Dikeman, M.Shifman and N.G.Uraltsev, Int. J. Mod. Phys. A 
12 (1997) 2075. 

[13] S.V.Mikhailov and A.V.Radyushkin, JETP Lett. 43 (1986) 712; Sov. J. Nucl. Phys. 
49 (1989) 494. 



12 



[14] M.Abramowitz and I.A.Stegun, Handbook of Mathematical Functions, New York 
(1965). 

[15] P.Colangelo, G.Nardulli and M.Pietroni, Phys. Rev. D 43 (1991) 3002. 

[16] PCca, P.Colangelo, G.Nardulli, G.Paiano and G.Preparata, Phys. Rev. D 26 (1982) 
1157. 

[17] PCea and G.Nardulli, Phys. Rev. D34 (1986) 1863. 

[18] M.Neubert, Phys. Rev. D 45 (1992) 2451; ibidem 46 (1992) 1076. 



13 



FIGURE CAPTIONS 

Fig. 1 

Plot of the ratio (O) as a function of E for several values of A (from top to bottom: 
A = 1.5,2.0,3.0 and 5.0 GeV). 
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